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Perturbations in the ionization fraction after recombination affect the Compton cooUng of density 
perturbations. Once the gas temperature starts to decouple from the CMB temperature, ionization 
fraction perturbations can have a significant infiuence on the subsequent gas temperature pertur- 
bation evolution. This directly affects the 21cm spin temperature of the gas, and also modifies the 
small-scale baryon perturbation evolution via the difference in baryon pressure. The effect on the 
gas temperature perturbations can be significant on all scales, and galactic-scale baryon perturba- 
tions are modified at the percent level at redshifts z > 100 where numerical simulations are typically 
started. 



I. INTRODUCTION 



Perturbations in the ionization fraction have almost no 
effect on the CMB power spectrum. This is because an 
isotropic photon field scattering from an arbitrary static 
anisotropic electron distribution will still look isotropic. 
Perturbed recombination only changes the CMB though 
effects such as the scattering of the perturbed temperature 
from the perturbed ionization fraction, which is a second- 
order process. However ionization fraction perturbations 
can influence the gas temperature perturbation evolu- 
tion at linear order by changing the local Compton cool- 
ing rate. This in turn will modify the spin tempera- 
ture perturbation relevant for dark-age 21cm observa- 
tions [U, 0, H, 0, § , and will also change the Jeans '-scale 
baryon density evolution due to the difference in baryon 
pressure. This can be significant at redshifts z > 100 
where numerical simulations are typically started, and 
should be included in any accurate calculation of the 
small-scale linear-theory baryon power spectrum at high 
redshift. 

The post-recombination effect of ionization fraction 
perturbations can be understood as follows. After the 
main recombination epoch the temperature is well below 
the ionization energy of the lowest states of hydrogen. 
The photoionization rate is therefore very low, and elec- 
trons captured by protons result in the net production 
of neutral hydrogen atoms. Recombination is therefore 
limited by the low probability for electron capture in the 
sparse gas, and the recombination rate depends predom- 
inantly on the atomic density and velocity distribution 
(temperature). Over-dense regions, which are also hot- 
ter due to adiabatic contraction, will have more recombi- 
nations, and hence have a lower ionization fraction than 
the background. This is turn reduces the coupling to 
the CMB temperature via Compton scattering because 
there are fewer free electrons. After recombination the 
gas temperature is below the CMB temperature due to 
adiabatic cooling; the residual lower ionization fraction 



in an overdensity therefore reduces the coupling to the 
higher CMB temperature, and hence makes the gas tem- 
perature lower than it would have been. A lower gas 
temperature implies a lower baryon pressure, and hence 
growth of small-scale baryon perturbations is less sup- 
pressed by pressure support. 

In this paper we give an approximate numerical analy- 
sis of perturbed recombination to quantify the linear ef- 
fect on the gas temperature and small-scale baryon den- 
sity evolution. We calculate the evolution only where 
linear-theory is expected to be quite accurate (redshift 
z > 30). Note however that the overall accuracy of the 
calculation is currently limited by the understanding of 
the complex background recombination process (see e.g. 
Refs. a detailed accurate treatment of perturbed 

recombination at all redshifts will have to wait until the 
detailed physics of recombination is more fully under- 
stood. 



II. PERTURBED RECOMBINATION 

The main effect of perturbed recombination during the 
dark ages can be obtained straightforwardly by consid- 
ering the ionization fraction evolution after the main re- 
combination event. This is because the perturbations 
grow after recombination, and the result at redshifts of 
a few hundred is not very sensitive to the initial value 
at the end of recombination. Defining an ionization frac- 
tion Xe, the number density of neutral hydrogen atoms is 
uhi = (1 — Xe)nH, where nn is the total number density 
of ionized and unionized hydrogen. We work with confor- 
mal time (derivatives denoted by an over-dot) and con- 
sider only scalar perturbations in the synchronous gauge 
(so the CDM peculiar velocity is zero). We assume a 
standard flat concordance cold dark matter cosmology 
in which the dark matter does not decay or annihilate 
(if it does the thermal history can be significantly al- 
tered [13]). 

At late times the evolution of the background ioniza- 
tion fraction is given by 
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where a is the scale factor and a is the recombination 
coefficient, a function of the gas temperature Tg fit by 



ac,{Tg/lQ''Kf 
l + c[Tg/lQ'^KY 



m^s 



(2) 



Here F is a fudge factor taken to be 1.14 from REC- 
FAST [ll| to caUbrate the result to a multi-level atom 
calculation, and — 4.309 x 10^^^, b = -0.6166, 
c = 0.6703, d = 0.5300. 

At late times the linear ionization fraction perturba- 
tion Ax then satisfies 
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where A^ are the fractional perturbations and 



A, 



b + c{Tg/ 10^ K)''{b-d) 
1 + c{Tg/lO^KY 



An 



(3) 



(4) 



Since the fractional baryon perturbation A„^ « A^ 
grows rapidly after recombination as the baryons fall in 
to the dark matter potential wells, the detailed evolution 
during recombination is not important for the late time 
amplitude of A^^ ■ The sign of A^^ is opposite to that of 
Ag because over-densities recombine more efficiently, so 
the ionization fraction is below that of the background 
at late times. 

Note that the overall calculation is limited by the pre- 
cision of the RECFAST model, in which a single fudge- 
factor accounts for deviations of an effective three level 
atom model from the full result. For example Ref. Q find 
a > 5% difference in at z ~ 200 compared to the REC- 
FAST model. Our results for the perturbation evolution 
are therefore of limited accuracy, though the relative im- 
portance of the ionization fraction perturbations should 
be roughly correct. 



III. GAS TEMPERATURE EVOLUTION 

After recombination the gas and CMB temperature re- 
main coupled for a while due to Compton scattering be- 
tween CMB photons and residual free electrons in the 
gas. However after some expansion the density and ion- 
ization fractions become lower and the gas starts to de- 
couple from the CMB and cool adiabatically. In detail, 
assuming purely Compton coolin g, t he background gas 
temperature evolves according to [111 [T^ 

+ 2nTg = ^7^7- (Tg - T,) , (5) 

3meC(l + UHc/nH + Xe) 

where ctt is the Thomson scattering cross-section, rifje 
is the number density of Helium, Tj is the CMB tem- 
perature, me is the electron mass and c is the speed of 
light. Gas temperature perturbations at wavenumber k 




FIG. 1: Evolution of the ionization fraction, matter and gas 
temperature perturbations with (solid) and without (dashed) 
accounting for ionization fraction perturbations. Perturba- 
tions are normalized to correspond to unit initial super- 
horizon curvature perturbation. The ionization fraction per- 
turbation can be positive (shown dash-dotted) during the 
main recombination event since overdensities can effectively 
delay recombination. The late-time residual local ionization 
fraction measured by A^^ is insensitive to the local timing of 
the main earlier recombination. 
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where we neglected helium fraction perturbations, v is 
the baryon velocity, and Ac is the fractional cold dark 
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matter perturbation. This equation depends explicitly on 
the ionization fraction perturbation A^.^ , which is often 
neglected. For example Ref. [l^ neglects the perturba- 
tions entirely, while Ref. [s*] gives an argument for them 
being much smaller than the density perturbations. This 
is true, but the effect on the temperature evolution can 
still be significant if not large. Ionization fraction per- 
turbations are modelled approximately in Ref. [10|] for 
the case of decaying and annihilating dark matter where 
the effects become large, but the impact on the standard 
CDM-scenario was not quantified. 

An overdensity has positive Aj-^ but recombines more 
fully than the background and hence has negative A^^j 
the additional effect of the ionization fraction perturba- 
tion is therefore to reduce the coupling to the CMB and 
hence slightly decrease the positive temperature pertur- 
bation in an overdensity. In addition the linear baryon 
density perturbations are slightly modified 



IV. BARYON EVOLUTION 

After recombination Compton scattering has no sig- 
nificant impact on the baryon density evolution, and the 
gas pressure is low so that terms involving p^/pg <C 1 can 
be neglected. The fractional synchronous gauge baryon 
density perturbation A^ then obeys the equation 
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where the sum is over the perturbations in the den- 
sity and pressure of the fluid components, just baryons 
and dark matter during matter domination. The baryon 
sound speed Cg is defined in terms of baryon density and 
pressure perturbations as = Spb/Spb- Assuming an 
ideal gas so that the gas pressure is given in terms of 
the density by Pg oc TgPg , the gas pressure perturbation 
4^9 = ^Pg/Pg is given by Q 
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X ([1 -I- Xe]{l - Ihc) + lHc"^p/mHc) , (8) 

where p is the mean particle mass and Ihc is the mass 
fraction in helium. This result must be used on scales 
where the baryon pressure is important [3]; the com- 
mon approximation that ~ dp/ dp is not accurate and 
also misses the linear effect of ionization fraction pertur- 
bations. The small-scale baryon perturbation growth is 
suppressed by baryon pressure, so the temperature per- 
turbation, which depends on the ionization fraction per- 
turbation, must be calculated properly for accurate re- 
sults on scales where kcg/'H ^ 1 just after recombination. 

Note that there may be important additional non- 
linear effects on the baryon evolution on ultra-small 



FIG. 2: Linear-theory transfer functions at redshifts z = 150 
(top plot) and a = 50 (bottom plot) for the density and tem- 
perature perturbations (top panel), and the fractional differ- 
ence due to inclusion of linear ionization fraction perturba- 
tions (bottom panel). Dashed lines show the result without 
including ionization fraction perturbations. Perturbations are 
normalized to correspond to unit initial super-horizon curva- 
ture perturbation. Note that non-linear effects will also be 
important on small scales at redshift z = 50. 



scales at recombination, e.g. see Refs. [11 [11] (though 
note that these papers incorrectly define = p/din 
the linear baryon evolution equation following Ref. [la]). 
In addition late-time non-linear evolution can be non- 
negligible on very small scales even at high redshift. 
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V. NUMERICAL RESULTS 

Numerical results for the gas temperature and baryon 
density can easily be obtained using perturbation evolu- 
tion equations implemented in standard Boltzmann codes 
such as CMBFAST [l^ and CAMB [H]. As of 2006 
neither code correctly uses Eq. ^ for the sound speed, 
largely because the error is harmless on currently ob- 
served scales. We have modified the CAMB code to 
do this correctly, including approximate ionization and 
gas temperature perturbation evolution equations as dis- 
cussed above^. For the ionization fraction perturbations 
we simply use a full set of perturbed RECFAST equa- 
tions (given in the Appendix of Ref. [5j]); these may not 
be accurate during recombination due to velocity escape- 
probability (and other) effects, but we have checked that 
the later-time dark-age ionization fraction perturbation 
is insensitive to these details. For redshifts z < 700 the 
result agrees well on most scales with integrating Eq. ^ 
from the start of recombination. Accuracy is likely to be 
limited by the accuracy of the background recombination 
model. 

Figure [T] shows the evolution of the various pertur- 
bations with time. For modes with wavenumber k = 
IMpc"^ the baryon pressure has almost no effect on the 
density perturbation evolution, but at k — 500Mpc~^ 
pressure suppresses the baryon perturbation growth. 
On both scales ionization fraction perturbations have a 
many-percent effect on the evolution of the gas tempera- 
ture fluctuations at redshifts z > 40 where the Compton 
coupling is not negligible. At 40 < z < 100 the fractional 
change in At^ is 5-10%. 

CoUisional coupling of the hydrogen 21cm spin tem- 
perature to the gas temperature means that changes in 
the gas temperature affect the dark-age 21cm fluctua- 
tions on all scales (for details see Ref. Q). At redshifts 
40 ^ z ^ 70 overdensities have a higher collisional cou- 
pling rate than the background and hence have a spin 
temperature closer to the gas temperature. Allowing for 
ionization fraction perturbations decreases the gas tem- 
perature in an overdensity and hence lowers the spin tem- 
perature further, giving a larger absorption signal. At 
z < 40 and the collisional spin-coupling rate falls rapidly 
with temperature; this can be seen as a change in the 
slope of the collisional de-excitation rate as a function of 
Tg for Tg < 30K (see e.g. Ref. (3). Here lowering the 
gas temperature due to ionization fraction perturbations 
leads to a lower collisional coupling rate, which actually 
causes the spin temperature to slightly increase, giving 
a smaller absorption signal than when A^,^ is neglected. 
The effect of ionization fraction perturbations perturba- 
tions is therefore to inc rease the signal at z > 40 and 

^ Code publicly available at |http : //camb Tlnf o/sources/ 1 



decrease it at z < 40. The 21cm power spectrum is mod- 
ified at the > 2%-level at z > 50, but the effect at z ~ 40 
is very small due to the cancellation of the competing 
processes. 

Figure O shows the perturbation transfer functions 
as a function of wavenumber. Ionization fraction per- 
turbations affect the gas temperature transfer function 
on all scales at the many-percent level. The effect on 
the baryons is small, but at the percent-level on scales 
where baryon pressure is important (i.e. modes inside 
the baryon sound horizon at recombination; for further 
discussion of baryon perturbation evolution see e.g. the 
Appendix of Ref. Q). Ionization fraction perturbations 
reduce the temperature in overdensities which reduces 
the pressure, giving a small increase in the baryon den- 
sity due to the smaller pressure support. The baryon 
density perturbations gravitationally influence the dark 
matter perturbations, though only at the O(10~^)-level, 
which is well below corrections due to non-linear evolu- 
tion on these scales. At later times the gas cools further 
and the baryon perturbations relax to follow the dark 
matter; the late-time linear-theory baryon perturbation 
becomes close to what it would have been without ion- 
ization fraction perturbations. 



VI. CONCLUSIONS 

Perturbation in the dark-age ionization fraction can 
have a significant effect on the gas temperature evolu- 
tion. This must be modelled for accurate calculation 
of the dark-age 21cm power spectrum. The effect on 
the 21cm power spectrum is at the 2%-level on all 
scales at redshift z ~ 50 (with a larger effect on very 
small scales), as discussed in more detail in Ref. [1]. 
Baryon evolution is also affected at the percent level 
on scales smaller than the baryon sound horizon at 
recombination, corresponding to wavenumbers larger 
than a few hundred inverse megaparsecs. Multi- fluid 
numerical simulations may need to use an accurate 
model of fluctuations in the gas temperature to achieve 
accurate results on galactic scales from a given set of 
cosmological parameters. 
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